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O ■ ABSTRACT 



We investigated the characteristics of the shallow decay phase in the early X- 

ray afterglows of GRBs observed by Swift X-Ray Telescope (XRT) during the 

Qh' period of January 2005 to December 2006. We found that the intrinsic break time 

O ■ at the shallow-to-normal decay transition in the X-ray light curve T^^^ is moder- 

^ • ately well correlated with the isotropic X-ray luminosity in the end of the shallow 

J^: decay phase (Lx.end) as T°,k = (9.39±0.64) x 10=^s(Lx,end/10^Vgs ■ s-i)-o-7i±o-03^ 

while T^j.]^ is weakly correlated with the isotropic gamma-ray energy of the prompt 

^ ■ emission E^^iso. Using T^^^^ — Lx,end relation we have determined the pseudo red- 

^D • shifts of 33 GRBs. We compared the pseudo redshifts of 11 GRBs with measured 

G> ■ redshifts and found the rms error to be 0.17 in log 2;. From this pseudo redshift, 

^^ we estimate that ~ 15% of the Swift GRBs have z > 5. The advantages of this 

distance indicator is that (1) it requires only X-ray afterglow data while other 

^ ■ methods such as Amati and Yonetoku correlations require the peak energy (Ep) 

of the prompt emission, (2) the redshift is uniquely determined without redshift 



% 



degeneracies unlike the Amati correlation, and (3) the redshift is estimated in 
advance of deep follow-ups so that possible high redshift GRBs might be selected 
for detailed observations. 
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1. Introduction 

The optical afterglow light curves of most GRBs show the smooth power-law decays 
with time (oc t~" with a typical index value a ~ 1). This is consistent with the prediction 
of the simplest model of a spherical blast wave propagating into a uniform medium, where 
the spectrum consists of several power-law segments in which F^ oc z/^'^t^" (Sari, Piran & 
Narayan 1998). 

However, after the advent of the HETE-2 satellite, prompt localization of GRBs sent 
to ground-based telescopes made it possible to observe the early afterglows which revealed 
that some GRBs show deviations from a smooth power-law light curve. For example, GRB 
021004 showed a highly variable light curve, with several bumps and wiggles from a simple 
power-law (e.g.. Fox et al. 2002; Mirabal et al. 2002). GRB 030329 also showed such a short 
timescale variability (Uemura et al. 2003; Sato et al. 2003). Furthermore, the early slow 
decay was observed in GRB 021004 and was interpreted by Fox et al. (2003) as arising from 
delayed shocks and continuous energy ejection from the central engine (Rees & Meszaros 
1998; Kumar & Piran 2000; Sari & Meszaros 2000). 

The early X-ray afterglow was found to be even more complex by the Swift observations. 
They generally consist of three distinct power-law segments (Nousek et al. 2006; Zhang et al. 
2006): (1) an initial steep decay with ai ~ 3, (2) a shallow decay with 02 ~ 0.5 and finally 
(3) a normal decay 0:3 ~ 1, where ai, 0:2 and 0:3 are power-low indices of temporal variability 
{a of t~"). The initial steep decay is most likely the transition from the prompt emission 
to the afterglow (Nousek et al. 2006; Zhang et al. 2006). Such a steep decay is due to the 
"curvature effect" of the high-latitude emission expected for the emission ceasing abruptly 
(e.g., Kumar & Panaitescu 2000; Dermer 2004; Dyks et al. 2005; Panaitescu et al. 2006; 
Yamazaki et al. 2006). In the normal decay phase, the data for many bursts are consistent 
with an ISM medium rather than a wind medium. However, for the shallow decay phase, its 
physical background has not been understood and it is the most mysterious feature in the 
early X-ray afterglows. 

As for shallow decay phase, Willingale et al. (2007) have studied the end-time of the 
shallow decay of the X-ray light curves of long GRBs, Ta. They suggested a possibility 
that Ta depends on the total energy of the outflow. Nava et al. (2007) have investigated the 
correlation between the end-time in the GRB frame T^ = Tg^/{l+z) and the isotropic gamma- 
ray energy of the prompt emission -E^^iso- They found that for the bursts in their sample, 
T° is weakly correlated with E^^iso- However, this correlation disappears when considering 
all bursts of known redshift and Ta. Liang et al. (2007) have also showed that there is no 
significant correlation between the break time between shallow and normal decay segments 
and the -E-y^iso- 
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In this paper we investigated the characteristics of the shallow decay phase in the early 
X-ray afterglows of GRBs observed by Swift X-Ray Telescope (XRT) during the period of 
January 2005 to December 2006. We found that the intrinsic break time at the shallow-to- 
normal decay transition in the X-ray light curve T^,.^ is correlated with the isotropic X-ray 
luminosity of the shallow decay phase, Lx- We tried to apply this relation to determine the 
redshifts of 33 GRBs and found that the distribution of the pseudo redshifts is similar to that 
of spectroscopically measured redshifts with more high pseudo redshift GRBs. Furthermore, 
we independently examine if there are any correlations among parameters of the prompt 
emission and the shallow decay phase for GRBs. We found that T°j.]^ is weakly correlated 
with the isotropic gamma-ray energy of the prompt emission -E^^iso- 

Throughout this paper, we adopt the cosmological parameters Hq = 71 km s~^ Mpc~^, 
Qu = 0.27 and ^a = 0.73. 



2. Data selection and analysis 

In this section, we present the data samples and the calculation method of parameters 
as Lx, E-y^iso and Tbrk- To examine the Lx and T°j,j^ relation, the measured redshift z and 
XRT data at the start and the end time (Tbrk) of the shallow decay phase are needed: The 
isotropic X-ray luminosity Lx is calculated from the X-ray spectrum and the break time at 
the shallow-to-normal decay transition Tbrk = (1 + ^)T^rk i^ obtained from the X-ray light 
curve. Here we define the shallow decay phase when the light curve decay index is flatter 
than the canonical value of a ~ 1. We found that 11 GRBs have well defined measured 
values of these parameters between January 2005 and December 2006. 

On the other hand, in order to examine the E^^iso — T^rk relation, the redshift z, the 
photon energy at the peak of the uF^, spectrum, Ep for calculating -E^,iso and Tbrk are needed. 
Ep is obtained from the GRB spectrum. A GRB spectrum is typically described by a Band 
function (Band et al. 1993), which is a smoothly-joint broken power- law characterized by two 
photon indices and Ep. We found that seven GRBs had measured value of these parameters 
between January 2005 and December 2006. Ep for these GRBs has been firmly identified 
from the Konus, HETE — 2 observations reported in Amati et al. (2006b). Generally it 
is difficult to determine Ep with Swift-BPCT data alone due to its narrow energy range of 
15 — 150 keV. For the GRBs whose Ep cannot be determined, we have developed a method 
to estimate Ep using the Ep — L^^iso correlation (Yonetoku et al. 2004; Ghirlanda et al. 2005) 
and additionally obtained Ep for 11 GRBs including GRB 050824 which has only the lower 
limit for Tbrk- 
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In this section we systematically analyze the XRT and the Burst Alert Telescope (BAT) 
data and show the details of calculation methods to obtain the parameters of Lx, E^^iso and 

^brk- 



2.1. XRT analysis 

2.1.1. Reduction 

The elapsed time of the events such as the break time is measured from the BAT trigger 
time in this paper. The XRT data presented here were obtained using the Window Timing 
(WT) and/or the Photon Counting (PC) modes (event grades 0—2 and 0—12 respectively). 
The XRT data were processed using xrtpipline into filtered event lists. Data were also 
filtered to eliminate time periods when the CCD temperature was warmer than — 50°C. 
These filtered data were then used to extract light curves and spectrum in the 0.5—10 keV 
energy range. For the PC data, the light curves and the spectrum are generally extracted 
from a circular region with a radius of 47" (the region size depends on Point Spread Function 
(PSF) of the X-ray afterglows.). The backgrounds are selected from an annulus region with 
radii of 94" and 188", excluding the X-ray source region near the GRB position. For the 
WT data, the light curves and the spectrum are extracted from a rectangular region of 94" 
by 47" . The background is selected from a rectangle of the same size as for the source region 
that is typically 47" away from the GRB position. In all cases, we used XSELECT to extract 
source and background and XSPEC version 11.3.2 to fit the spectra. 

The data obtained in the PC mode sometimes suffered from pile-up when the observed 
sources were brighter than 0.5 cts/s. These data were corrected for pile-up by adopting the 
method described in Vaughan et al. (2006). We used an annular extraction region, with a 
9.4" inner radius and a 47" outer radius. In order to determine the correction factor for the 
annular aperture, we modeled the PSF using XIMAGE. The effective area is corrected using 
the calibration data and xrtmkarf. 



2.1.2. Estimation ofTy^rk 

We consider 21 GRBs with known redshift and identified the shallow decay phase in the 
X-ray light curves. In order to determine Tbrk, we fitted the X-ray light curve obtained in 
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Section 2.1.1 to a broken power-law model: 



No X {t - toy' {t<T, 



Nit) 

-to 



brkJ 



No X (Tbrk - to)"'' X ji^^ ) (t > Tbrk), 



-a2 



where No is the normalization in units of counts s~^, ai and a2 are the temporal power- 
law indices before and after the break time Tbrk- The temporal decay slopes and the break 
times of the shallow-to-normal transition are summarized in Table 1. These parameters were 
determined using the BAT trigger time as to. 



2.1.3. Estimation of Lx 

In order to determine Lx, we fitted the XRT spectra obtained in Section 2.1.1 with the 
single power-law {N{E) ex E~^, where F is the differential photon index). We here adopt a 
k-correction to Lx- Let us write N{E)dE = AE~^dE, then Lx integrated from E^ to E^ in 
the GRB frame is given by 

i'po\2-r _ /po\2-r 
Lx = 4nd,{z)\l + .)"'M ^^"^ ^^_[^'^^ , (2) 

where dj, is the luminosity distance, F is the power-low index, E^ and E^ are measured in 
the GRB frame. To determine E^ and E^ we consider GRB at the Swift average of z ~ 3 
observed in 0.5 — 10 keV band so that E^ = 2 keV and E^ = 40 keV. We tested for the 
X-ray luminosity at initial {Lx,im), median (Lx,med) and end {Lx,end) of the shallow decay 
light curve. The results are summarized in Table 2. Some events did not have observation at 
beginning of the shallow decay phase for convenience of Swift observation. For these events, 
we show only the parameters at the end of the shallow decay light curve. 



2.2. BAT analysis 

2.2.1. Reduction 

We analyze the BAT data of the Swift GRBs observed during the period January 2005 
to December 2006. The BAT data for the GRB samples were processed using standard 
Swift-BAT analysis software as described in the BAT Ground Analysis Software Manual 
{http : / /heasar eg sfc.nasa.gov /docs /swift /analysis/). Each BAT event was mask-tagged 
using batmaskwtevt task with the best fit source position. All of the BAT spectra have been 
background subtracted with this method. 
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We extracted the spectra in the energy range 15 — 150 keV over the period corresponding 
to T90 excluded during slew. If the satellite started to slew during T90, we extracted the 
spectrum before the slew start times. If the satellite did not slew during T90, we extracted 
the spectrum over the whole Tqq. All spectra were fitted with XSPEC version 11.3.2. The 
detector responses are generated from batmaskevt. 



2.2.2. EsUmaUon of E, 



7, ISO 



In this section, we consider 11 Swift GRBs for which spectroscopic redshifts and Tbrk 
are both available. A GRB spectrum is typically described by a Band function (Band et al. 
1993). The typical values of two photon indices are Fi ~ 1.0 and r2 ~ 2.2, respectively. 
However, the spectral peak energy values of GRB is typically E^ ~ 250 keV, i.e., above the 
BAT energy band (15—150 keV). Actually most spectra observed by the BAT are well fitted 
by a single power-law function. Generally, with the BAT observations alone, we cannot 
determine E^ and the high energy photon index r2. 

We estimate E^ by using the E-^ — L^^iso relation (Yonetoku et al. 2004; Ghirlanda et al. 
2005) where L^^iso is the gamma-ray isotropic luminosity ^. Since the photon index obtained 
from the BAT is distributed < 2 in most case, the peak energy Ep of bursts is expected to 
be above the BAT energy band. Therefore, we assume a broken power-law shape for the 
spectra time-averaged over the GRB duration: 

^^^^ -\Ax Ep^'^-'^^E-^^ for E > E,. ^^^ 

where A is a normalization, Fi and F2 are the low and high energy photon indices, respec- 
tively, and Ep is the peak energy. A and Fi can be determined from BAT data, and F2 is 
assumed to have the typical value of 2.2. The time-averaged flux in a given bandpass {Ei, 
E2), where Ei and E2 are the minimum and the maximum energy, respectively, as a function 
of Ep can be given as: 

Fe,-e,{Ep)= / EN{E)dE. (4) 

J El 

Therefore, the isotropic luminosity L-y,iso is calculated as 

L^,iso = Airdlx Fei-E2{Ep), (5) 



^The Ep — -E-y^iso Amati relation (Amati et al. 2002; Friedman & Bloom 2005; Amati 2006) is a similar 
luminosity relation. However, this does not produce a unique or necessarily well-determined Ep value because 
there is an intrinsic redshift degeneracy in the Amati relation (Li 2006; Schacfcr & CoUazzi 2007). 
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where di, is the luminosity distance. 

Meanwhile, the Ep — L^^iso correlation was proposed by Yonetoku et al (2004). Ghirlanda 
et al. (2005) re-examined this correlation with an enlarged sample and showed the correla- 
tions with smaller scatter: 

P / r ■ \ 0.48±0.01 



(4.88±0.06)xf _-^7;- ) . (7) 



100 keV ' ■ ' ' V1.9 X W^ J 

From equations (6) and (7), we can obtain the peak energy E^. Then we apply the 
-Ep — -E-y,iso relation (Amati 2006b) to estimate E^^[so, by assuming Ep calculated above. 
For the bursts with observed -E^,iso values, we compared -E^,iso calculated from our method 
and observed value. Then we confirmed that calculated values are approximately consistent 
with the observed values. In Figure 1, the dotted lines show difference with a factor of two 
between the observed and calculated values. We see that except for short GRB, the observed 
and calculated values agree within a factor of two. We fitted the BAT spectra obtained in 
Section 2.1 with the single power-law {N{E) oc E~^, where F is the differential photon 
index). Then we apply the Ep calculation method to these BAT photon indices. The results 
are summarized in Table 3. 



3. Results 

3.1. Lx — ^brk relation 

Figure 2 shows the distribution of the intrinsic break time at the shallow-to-normal 
decay transition in the X-ray light curve in the GRB frame T^^.^ as a function of the X- 
ray afterglow luminosities at different epochs: Lx,im at the beginning of the shallow decay, 
Lx,mcd at the median epoch and Lx,cnd at the end. Eleven GRBs in our sample have redshift 
measurements and XRT data in the start and the end time (Tbrk) of the shallow decay phase. 
There are moderately good correlations between afterglow luminosities and Tbrk- In these 
samples GRB 060607A is unusual since it has an abrupt break at Tbrk similar to GRB070110 
(Troja et al. 2007) and the origin of shallow decay could be different from others (Liang et 
al. 2007). Considering GRB060607A as an outlier, the correlation coefficient are 0.47, 0.65, 
0.65 for 9 degrees of freedom (10 GRBs) for Figure 2 a, b and c, respectively. The chance 
probability are 0.17, 0.042, 0.042, respectively. 

When we adopted the power-law model to the Lx,mcd - T^rk ^^^ ^^^ Lx,end - T^rk relation. 



the best-fit function are 



J- K, 



Tbrk , , _ (10-6 ± 0.78) X 10^ (j^^r^^-T ) for Lx,^ed, 



-0.68±0.03 



brk — /-I , \ l^; — \ / r \ -0-71±0.03 K^ ) 

^^^^^ [(9.39±0.64)xl03(^^p4^^j for Lx,end. 

Figure 3 shows Lx,cnd — ^brk relation for all the bursts that have Tbrk measurements 
including those missing the observational data at the beginning of the shallow decay phase. 
We found that two types of the bursts seem to deviate from the Lx,end — ^brk correlation: 
the bursts which have small X-ray luminosity at Tbrk, and the bursts which have abrupt or 
chromatic X-ray light curve breaks (GRB 060607A, 050319, 050401 (Panaitescu et al. 2006), 
060210 (Stanek et al. 2007), and 060927 (Ruiz-Velasco et al. 2007)). These plots are shown 
as open circles. 

In Figure 4, we show the integrated energy in the shallow phase {Eg/E^^iso) and E^^iso- 
We see that Eg is typically (0.01~0.1) -E^.iso- This result is similar to Figure 3-d of Liang et 
al. (2007). 



3.2. -E-y^iso — Tbrk relation 

Figure 5 shows the distribution of T°j,j^ and the isotropic gamma-ray energy of the prompt 
emission E^^iso- We found that T^^^. is weakly anti-correlated with -E^,iso in logarithmic scale 
^. The correlation coefficient is 0.49 for 16 degrees of freedom (15 GRBs except GRB 050824 
which has only the lower limit for T^j.]J and the chance probability show value of 0.057. This 
correlation suggests that the larger the isotropic equivalent energy the earlier the end time 
of the shallow decay phase. 

Note that Nava et al. (2007) suggested that Ta weakly correlates with -E^^iso in their 
sample. However, this correlation disappears when considering all bursts of known redshift 
and Ta. Ta is obtained by fitting the X-ray light curves with the prompt emission and the 
afterglow component functions (Willingale et al. 2007). In their methods, even if the shallow 
decay phase is not clearly observed in the X-ray light curve, they can obtain the Ta value. 
On the other hand, we defined the shallow decay phase when the light curve decay index 
is flatter than the canonical value of a ~ 1. As for the GRBs satisfying this condition, our 



^In order to check validity of that assumption of r2 = 2.2, we tested the case of r2 =2.1 or 2.5. Ahhough 
the slope of the plot changed with the r2 value, it did not affect the correlation between S^^so and T^^^^. 
We also plotted only the data with firmly identified i?-y,iso- Though the dispersion of the data point is large, 
there seem to be a trend that the larger i?-y,iso, the earlier T^^^^. 
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plot and the plot showed in Nava et al. (2007) are approximately consistent. Furthermore, 
Liang et al. (2007) showed that there is no significant correlations between T^^.^^ and -E^,iso 
in their sample. They defined the shallow decay phase has a < 0.7 since the decay slope 
of the normal decay phase predicted by the external GRB models is generally greater than 
0.7. Therefore the GRBs used for the T^^.^^ — -E^,iso correlation test are slightly different from 
ours. 



4. Estimation of the redshift using Lx,cnd — ^brk correlation 

In the previous section we found a moderately good correlation between T^^y. and the 
isotropic X-ray luminosity in the end of the shallow decay phase {Lx,cnd) as well as the 
weak correlation between T^^,]^ and E^^iso- The correlation coefficients are not so good to 
insist the relations. Nevertheless, in this section assuming that Lx,end — T'brk correlation is 
correct we apply the relation to determine the redshifts of GRBs observed by Swift X-Ray 
Telescope (XRT) during the period of January 2005 to December 2006. This is challenging 
and important since in ~ 200 GRBs observed by Swift only ~ 50 have the spectroscopically 
measured redshifts so that another distance indicator using only BAT and XRT data is of 
great value. 

We first rewrite equation (8) using only the observed quantities and the unknown redshift 

z as 

(1 + zr^-^^hizY-^' = 4.78 X 10-^ X T,-^ X [a x ^^"^ _ [.^'^^ ) , (9) 



where 



r{z) = r , "^^ ^=. (10) 

Jo ^Qmi^ + ZY + Q^ 

Since the left hand side of equation (10) is a monotonically increasing function from 
zero, there is only one solution of z for any observed values of Tbrk and -Fx.cnd- We call the 
redshift obtained by this method as the pseudo z from Lx,cnd — ^brk correlation. 

In Table 4 we show the list of the pseudo redshifts obtained for GRBs observed by 
Swift X-Ray Telescope (XRT) during the period of January 2005 to December 2006. Figure 
6 compares the pseudo redshifts with spectroscopically determined redshifts. Although the 
error bars of the pseudo redshifts are rather large, we see that the pseudo redshift determined 
by our method is a relatively good measure of the redshift of GRB for which no spectroscopic 
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information is available. The correlation coefficient is 0.58 for 15 degrees of freedom (16 
GRBs). 

In Figure 6, the dotted lines show difference by a factor of two between the observed 
and pseudo redshifts. We see that except for two GRBs, which have abrupt/chromatic X-ray 
light curve breaks (GRB 060210 and 060607A: different by a factor of three), the observed 
and pseudo redshifts agree within a factor of two so that we may say that our pseudo redshift 
is the measure of the redshift if we allow a factor of two error. 

In Figure 7 we show the cumulative distribution of the spectroscopically measured red- 
shifts ((1), the dashed line), the pseudo redshifts for GRBs with no spectroscopically mea- 
sured redshifts ((2), the solid line) and the total ((l) + (2), the dotted line). Note that the 
normalization of the cumulative distribution is the total one for (1) and (2). The mean red- 
shifts are 2.2 and 2.6 for the observed and the total redshifts, respectively. There is a slight 
difference in the distribution between of the pseudo z and observed z, but it is probably 
due to the difficulty of obtaining spectra at high redshifts. Figure 7 suggests that ~ 15% 
of GRBs have redshifts greater than 5. This is consistent with the constraints from the 
optically observed GRBs (Tanvir & Jakobsson 2007) 



5. Discussions 

Since T°j,j^ — -E^^iso correlation is weak we here argue mainly the physical implications of 
T^j,^ — Lx,cnd correlation as described by T^^.^. oc -^^^end- Here we note it is reasonable that 
T^^y. — Lx,cnd correlation is better than T^^^.^^ — Lx,im correlation since in the shallow decay 
phase more energy is emitted in the end phase than the initial phase. Interestingly this 
correlation is similar to the burning time tn and the luminosity L relation of the hydrogen 
main sequence star {tn oc L~°'^). From the theory of stellar evolution L and tn mainly 
depend on the mass of the star M as 

L(xM\ (11) 

tH oc — oc M'^. (12) 



Eliminating M from above two equations we have tn oc L 



-0.67 



5.1. Energy Injection Model 

The energy injection model (Nousek et al. 2006; Zhang et al. 2006; Granot & Kumar 
2006) is thought to be that energy is injected continuously into the external shock so that 
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the flux decay becomes slower than the usual oc t^^. The injection may be caused by (1) a 
long-lived central engine or (2) a short-lived central engine ejecting shells with some range 
of Lorentz factors. 

First we consider (1), the long-lived central engine model. This scenario requires the 
central engine to remain active until the end of the shallow decay phase [T^j.]^), which is in 
many cases 1,000 — 10,000 s. Since the X-ray afterglow is a good indicator of the kinetic 
energy, the central engine injects energy with the kinetic luminosity proportional to the X- 
ray luminosity L^in oc Lx oc t~"^ The kinetic energy (-Ekin) of the afterglow increases as a 
function of time oc t^~"\ and the observed correlation suggests that the total kinetic energy 
is anti-correlated with the lifetime of the central engine, 

^kin,c„d OC Lx,onXri. OC (T^^) "^+' ~ (T^V) ""•'' . (13) 

where -Ekin.end is the kinetic energy of afterglow in the end of the shallow decay phase. 

Next we consider (2) a short-lived central engine with some range of Lorentz factors 
of ejected shells. After the internal shocks, shells are rearranged so that outer shells are 
faster and inner shells are slower. This configuration may also occur if the central engine 
ejects faster shell earlier. Outer shells are slowed down by making the external shock. Once 
the Lorentz factor of the leading shocked shell drops below that of a following slower shell, 
the slower shell catches up with the shocked shell, injecting energy into the forward shock. 
Since the Lorentz factor of the afterglow is proportional to F oc E^^^T~^^^n~^/^, the observed 
correlation suggests that the Lorentz factor of the slowest shell, which has almost all energy, 
is nearly proportional to the total kinetic energy, 

Fslow (X Ei°,end^"'/', (14) 

where n is the ambient density. 



5.2. Inhomogeneous Jet Model 

In the inhomogeneous jet model, it is assumed that we observe more energetic com- 
ponents in the GRB jet at later times as the external shock decelerates and the visible 
region increases. The shallow decay phase is produced by the superposition of the afterglow 
emission from the off-axis components. This phase ends when the whole jet is observed in 
the ring-shaped jet model (Eichler & Granot 2006) or when the mini-jets merge and the 
inhomogeneities are averaged out in the multiple mini-jets model (Toma et al. 2006). 

In the ring-shaped jet model, the X-ray luminosity at the end of the shallow phase is 
given by Lx,cnd oc -Ekm,iso(^brk)~^' ^^^ the opening angle of the whole jet is determined by 
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9j oc -E'kin iso(^brk)^^^^^''^) where -Ekin.iso is the isotropic kinetic energy of the whole jet. Thus 
the observed correlation suggests 9j oc -^kinfso^^^^- Since the collimation-corrected kinetic 
energy is calculated as -Ekin oc -Ekin,iso^?! we obtain i^kin oc Om^^^. 



5.3. Time-dependent Microphysics Model 

This model considers that the microphysical parameters, such as the energy fraction 
that is shared to electrons ec and magnetic field ee, depend on time (loka et al. 2006; Fan & 
Piran 2006; Panaitescu et al. 2006b). We usually assume that the micro-physical parameters 
do not vary and in fact, constant eg and cb are consistent with the observation of late time 
afterglows (Yost et al. 2003). However, the behavior of these parameters in the early time 
afterglow is not yet known. 

In the model, the microphysical parameters vary in the early afterglow. After reaching 
the equipartition value, the microphysical parameters remain constant as observed in the 
late time afterglow. The X-ray luminosity Lx is given by the bolometric kinetic luminosity 
L as Lx ~ GcL. Since L oc -Ekin^"^, the shallow X-ray light curve Lx oc t~^/^ suggests that €e 
evolves as ec oc t^'^ (loka et al. 2005). The observed correlation suggests that the saturation 
occurs at the Lorentz factor Fsat oc -^kin cnd^^^^^ ^^ in equation (14). 



5.4. Pulsar Model 

Troja et al. (2007) showed that the abrupt drop of the X-ray light curve observed in 
GRB 070110 cannot be explained by an external shock as the origin of the shallow decay 
phase and implies the long-lived central engine. Furthermore, they suggest that the shallow 
decay phase might be powered by a spinning down central engine, possibly a millisecond 
pulsar. Motivated by this suggestion let us consider that the slowly changing usual shallow 
decay is also due to the pulsar activity. If the dipole magnetic field is constant, the luminosity 
of the pulsar decreases as oc t^^, so that we need the increase of the dipole magnetic field to 
interpret the shallow decay phase. From the total energy conservation, we have 

d /I . B^R^\ 

where we assume that the rotational energy goes into the magnetic energy by some unknown 
mechanism satisfying the observed luminosity. 
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Integrating the above equations from the beginning of the shallow decay phase to to T^j.^, we 
have 



l^'n2 n2, , {B'-Bi)R^ _ ^ . 1 



-/(fi - r^o) H - --^x,cnd^brkY" 



1 



1 \ 1-ai 



-'brk 



;i7) 



In this model T^^.^^ is essentially determined by the total energy conservation as 

hnl ~ Lx,endT°,k/(l - «l) = i^s. (18) 

Therefore the energy of the shallow decay phase Eg is essentially the total energy of the 
initial pulsar. Figure 4 suggests that the initial rotational period is lms~ 3ms, which is an 
appropriate value in this model. Using equation (16) and T°j,i^-Lx,end relation we have 

BK^[^y\TLr^- (19) 

Since Qq is known from Eg, the above equation tells us that the initial strength of the 
magnetic field determines T^^^^. 



6. Summary 

From our observational results, we found that the intrinsic break time at the shallow- 
to-normal decay transition in the X-ray light curve T^^.^^ is moderately well correlated with 
the isotropic X-ray luminosity in the end of the shallow decay phase (Lx,cnd) as T°j,j^ = 
(9.39 ± 0.64) X 103s(Lx,cnd/10^^ergs ■ s-i)-o-7i±o-03, while T^V is weakly correlated with the 
isotropic gamma-ray energy of the prompt emission -E^^iso- Using this relation we have 
determined the pseudo redshifts of 33 GRBs and found that the distribution of the pseudo 
redshifts is similar to that of spectroscopically determined redshifts. Since the T^^y- — Lx,end 
relation does not have an intrinsic redshift degeneracy, we can determine the redshift of the 
GRB uniquely. The T^j.]^ — Lx,end relation does not require the parameters of the prompt 
emission so that it may be useful to determine the redshift of Swift GRBs since the energy 
band of Swift is typically below the peak energy of the prompt emission. Our results suggest 
that ~ 15% of GRBs have z > 5. This means an exciting possibility such that the redshift 
is estimated in advance of deep follow-ups and possible high redshift GRBs {z > 6.3) might 
be selected for detailed observations and identified finally in near future . 

We discussed the implications of the T^^.^^ — Lx,end relation for some theoretical models 
recently proposed to explain the shallow decay light curve. In each model, we obtain an 
additional condition for the models to be satisfied from the T^^.^^ — Lx,cnd relation. Other 
models including two-component jet (Granot & Kumar 2006; Jin et al. 2007), dust scattering 
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(Shao & Dai 2007), and relativistic wind bubbles produced by the interaction of an ultra- 
relativistic electron-positron-pair wind with an outwardly expanding fireball (Dai 2004; Yu & 
Dai 2007), have also been proposed, but the detailed discussion for these models are beyond 
the scope of this paper. 

We are grateful to T. Sakamoto, G. Sato, and all the members of the Swift Team for 
their technical guidance of the analysis of Swift BAT and XRT data. R.S. is supported 
by the Research Fellowships for Young Scientists (2005-2007) of the Japan Society for the 
Promotion of Science. This work is supported in part by the Grant-in-Aid from the Ministry 
of Education, Culture, Sports, Science and Technology (MEXT) of Japan, No. 18740147 
(K.I.) andNo.l9540283,No.l9047004,No. 19035006(T.N.). 



REFERENCES 

Amati, L., et al. 2002, A&A, 390, 81 

Amati, L. 2006, MNRAS, 372, 233 

Amati, L. 2006b, submitted to II Nuovo Cimento C (astro-ph/0611189) 

Band, D. L., et al. 1993, ApJ, 413, 281 

Dai, Z. G. 2004, ApJ, 606, 1000 

Dermer, C. D. 2004, ApJ, 614, 284 

Dyks, J., Zhang, B., & Fan, Y. Z. 2005, ApJsubmitted (astro-ph/0511699) 

Eichler, D., & Granot, J. 2006, ApJ, 641, L5 

Fan, Y. & Piran, T. 2006, MNRAS, 369, 197 

Fox, D. W., Blake, C., & Price, P. 2002, GCN Circ. 1470 

Fox, D. W., et al. 2003, Nature, 422, 284 

Friedman, A. S., & Bloom, J. S. 2005, ApJ, 627, 1 

Ghirlanda, G., Ghisellini, G., Firmani, C., Celotti, A., & Bosnjak, Z. 2005, MNRAS, 360, 
L45 

Granot, J., & Kumar, P. 2006, MNRAS, 366, L13 



-15- 

loka, K., Kobayashi, S., & Zhang, B. 2005, ApJ, 631, 429 

loka, K., Tonia, K., Yamazaki, R., & Nakamura, T. 2006, A&A, 458, 7 

Jin, Z. P., Yan, T., Fan, Y. Z., & Wei, D. M. 2007, ApJ, 656, L57 

Kumar, P., & Piran, T. 2000, ApJ, 532, 286 

Kumar, P., & Panaitescu, A. 2000, ApJ, 541, L51 

Mirabal, N., Halpern, P., Chornock, R., & Filippenko, A. V. 2002, GCN Circ. 1618 

Li, L.-X. 2006, MNRAS, 374, L20 

Liang, E.-W., Zhang, B.-B., & Zhang, B. 2007, ApJsubmitted (astro-ph/07051373) 

Nava, L., Ghiselhni, G., Ghirlanda, G., Cabrera, J. L, Firmani, C, & Avila-Reese, V. 2007, 
MNRAS, 377, 1464 

Nousek, J. A., et al. 2006, ApJ, 642, 389 

Panaitescu, A., Meszaros, P., Gehrels, N., Burrows, D., & Nousek, J. 2006, MNRAS, 366, 
1357 

Panaitescu, A., Meszaros, P., Burrows, D., Nousek, J., Gehrels, N., O'Brien, P., & Wilhngale, 
R. 2006b, MNRAS, 369, 2059 

Rees, M. J., & Meszaros, P. 1998, ApJ, 496, LI 

Ruiz-Velasco, A. E. 2007, ApJsubmitted (astro-ph/07061257) 

Sari, R., Piran, T. & Narayan, R. 1998, ApJ, 497, 17 

Sari, R., & Meszaros, P. 2000, ApJ, 535, L33 

Sato, R. Kawai, N., Suzuki, M., Yatsu, Y., Kataoka, J., Takagi, R., Yanagisawa, K., & 
Yamaoka, H. 2003, ApJ, 599, L9 

Schaefer, B. E., & Collazzi, A. C. 2007, ApJ, 656, L53 

Shao, L., & Dai, Z. G. 2007, ApJ, 660, 1319 

Stanek, K. Z., et al. 2007, ApJ, 654, 21 

Tanvir, N. R., & Jakobsson, P. 2007, Phil. Trans. Roy. Soc. A (astro-ph/0701777) 



-16- 

Toma, K., loka, K., Yamazaki, R., & Nakamura, T. 2006, ApJ, 640, L139 

Troja, E., et al. 2007, ApJ, 665, 599 

Uemura, M., et al. 2003, Nature, 423, 843 

Vaughan, S., et al. 2006, ApJ, 638, 920 

Willingale, R., et al. 2007, ApJ, 662, 1093 

Yamazaki, R., Toma, K., loka, K., & Nakamura, T. 2006, MNRAS, 369, 311 

Yonetoku, D., Murakami, T., Nakamura, T., Yamazaki, R., Inoue, A. K., & loka, K. 2004, 
ApJ, 609, 935 

Yost, S. A., Harrison, F. A., Sari, R., & Frail, D. A. 2003, ApJ, 597, 459 

Yu, Y. W., & Dai, Z. G. 2007, A&A, 470, 119 

Zhang, B., Fan, Y. Z., Dyks, J., Kobayashi, S., Meszaros, P., Burrows, D. N., Nousek, J. A., 
& Gehrels, N. 2006, ApJ, 642, 354 



This preprint was prepared with the AAS I^TJtjX macros v5.2. 



-17- 



T3 

IP 



3 







• long GRB 







./\/\ 


o 
o 

o 

o 

o 




° short GRB 


-0 


^ 




o 


r 


- 




:.^ 




\ 


o 
o 

•0 


r 




/ ■' 






-, 


o 




/\ 








■ 



-10 ' 10 " 0.01 0.1 1 

Observed E. 



10 



100 1000 



10' 



JlO^"^ ergs) 



Fig. 1. — Comparison of the isotropic gamma-ray energy of the prompt emission -E^^iso 

see Section 2.2.2). 
The dotted hnes show difference with a 
factor of two between between the observed and calculated values. 



derived from our calculation method with firm estimated values -E°^so 
The dashed line is the calculated -E^,iso = -^'"liso- 
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Fig. 2. — Distribution of the intrinsic break time at the shallow-to-normal decay transition 
in the X-ray hght curve T^^.^. = Tbrk/(1 + 2;) as a function of the isotropic X-ray luminosity 
of the shallow decay phase Lx at different epochs: (a) Lx,im at the beginning of the shallow 
decay, (b) Lx,med at the median epoch and (c) Lx,end at the end. The dashed line is the 
best-fit power-low model for the data (Eq. 8). The open circle shows the unusual afterglow 
of GRB 060607A. 



-19- 




L , of shallow decay phase (10 erg/s) 



Fig. 3. — Distribution of the intrinsic break time at the shallow-to-normal decay transition 
in the X-ray hght curve T^^.^. = Tbrk/(1 + z) as a function of the isotropic X-ray luminosity 
at the end of the shallow decay Lx,cnd ■ The open circles show the unusual afterglows which 
have an abrupt break at Tbrk or a chromatic X-ray light curve break (see Section 3.1). 
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Fig. 4. — Distribution of the integrated energy in the shallow phase [Eg) and the isotropic 
gamma-ray energy of the prompt emission E^^is^- Eg is typically (O.01~0.1)-E.^^iso- 
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prompt emission -E'.^^iso- We found that T°j.j^ is weakly anti-correlated with -E^^^so in logarithmic 
scale. 
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Fig. 6. — Comparison of pseudo redshifts and the observed values except for the burst which 
have small X-ray luminosity at Tbrk- The pseudo redshifts were estimated from the correlation 
between the isotropic X-ray luminosity at the end of the shallow decay Lx,enA and the X-ray 
break times at the shallow-to- normal decay transition in the GRB frame T^^.-^- The open 
circles show the unusual afterglows, which have abrupt/chromatic X-ray light curve breaks. 
The dashed line is the pseudo redshift = measured redshift. The dotted lines show factor of 
two difference between the observed and the pseudo values. 
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Table 1: Temporal Parameters of the shallow-to- normal decay phase in the X-ray light curves 
for known-redshift GRBs. Parameters of tgta and tsto show the fitting range and the fitting 
model is a broken power-law (Eq. 1). 



GRB 



t 



sta 



^sto 

(s) 



ai 



logTbrk 

(s) 



a2 



Reduced x^ 
(d.o.f) 



050319 


3.24 


389 


417809 


0.49±0.03 


4.36±0.10 


1.27±0.26 


1.02 


;2i) 


050401 


2.90 


1020 


20108 


0.43±0.08 


3.73±2.39 


1.42±0.55 


1.01 


;i2) 


050505 


4.27 


2832 


44835 


0.15±0.15 


3.88±2.90 


1.17±0.05 


0.62 


;i7) 


050814 


5.3 


5759 


215431 


0.63±0.04 


4.89±0.06 


1.99±0.82 


1.19 


;i5) 


050824 


0.83 


6091 


603407 


0.52±0.04 


- 


- 


1.99 


(9) 


051016B 


0.9364 


4055 


278391 


0.73±0.06 


4.65±0.10 


1.36±0.16 


0.68 


;22) 


060115 


3.53 


749 


209997 


0.70±0.04 


4.73±0.10 


1.80±0.53 


1.10 


;24) 


060202 


0.783 


4942 


745686 


0.87±0.02 


- 


- 


1.07 


;46) 


060210 


3.91 


3849 


896348 


0.92±0.07 


4.47±1.82 


1.37±0.03 


0.97 


;38) 


060502 


1.51 


281 


676831 


0.54±0.05 


4.39±3.71 


1.19±0.05 


0.99 


;i7) 


060526 


3.221 


836 


220231 


0.41±0.06 


4.18±3.75 


1.46±0.24 


0.77 


;i3) 


060604 


2.68 


3470 


208383 


0.35±0.08 


4.31±3.38 


1.28±0.09 


0.99 


;i8) 


060605 


3.78 


230 


74165 


0.42±0.04 


3.93±2.54 


2.01±0.07 


1.14 


;25) 


060607 


3.08 


493 


48519 


0.42±0.02 


4.11±2.28 


3.62±0.07 


2.12 


;36) 


060614 


0.125 


4430 


555793 


0.11±0.04 


4.85±3.70 


2.24±0.13 


2.20 


;36) 


060714 


2.71 


283 


284999 


0.47±0.15 


3.59±3.10 


1.20±0.04 


0.98 


;27) 


060729 


0.54 


681 


180616 


0.09±0.01 


4.72±3.28 


1.18±0.05 


- 




060906 


3.686 


404 


139753 


0.13±0.08 


4.01±0.04 


1.76±0.13 


0.74 


;i2) 


060908 


2.43 


80 


92436 


0.67±0.05 


2.81±0.01 


1.41±0.04 


1.08 


;27) 


060927 


5.6 


70 


191186 


0.63±0.20 


3.51±3.25 


1.69±0.31 


0.57 


;io) 


061121 


1.314 


192 


92548 


0.30±0.23 


3.51±2.27 


1.25±0.02 


1.35 


;67) 
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Table 2: GRB samples used to study the X-ray luminosity and energy in the shallow decay 
phase. Lx is given in units of 10^^ erg s~^, and Eg is given in units of 10^^ erg s~^. 



GRB 


r 


Lx,im 


-^X.mcd 


-^X.cnd 


^s 






{W^ erg/s) 


(10^^ erg/s) 


(10^^ erg/s) 


(IQSi erg/s) 


050319 


2.21±0.15 


64.7±15.3 


9.82±2.41 


7.24±1.77 


6.69±1.79 


050814 


2.17±0.15 


27.5±6.3 


2.55±0.52 


1.67±0.32 


4.45±1.06 


060115 


1.84±0.09 


10.4±2.7 


0.865±0.251 


0.539±0.156 


1.75±0.49 


060502 


2.07±0.13 


7.77±2.79 


0.922±0.215 


0.639±0.155 


1.20±0.45 


060526 


2.35±0.17 


10.1±2.8 


4.11±0.91 


3.17±0.78 


1.63±0.64 


060605 


2.29±0.13 


55.1±15.5 


16.6±2.7 


12.5±2.1 


3.30±0.96 


060607A 


1.80±0.09 


186.±30. 


62.4±8.9 


47.3±6.2 


21.8±3.5 


060714 


2.91±0.17 


27.3±7.0 


11.0±2.0 


8.28±1.46 


1.18±0.41 


060729 


2.16±0.07 


0.223±0.049 


0.160±0.024 


0.150±0.022 


0.554±0.122 


060906 


2.29±0.22 


8.55±3.35 


6.10±1.54 


5.61±1.35 


1.31±0.52 


061121 


2.21±0.08 


31.3±5.3 


16.8±2.6 


14.0±2.3 


2.29±0.40 


050401 


2.10±0.06 


- 


- 


61.8±8.4 


- 


050505 


2.09±0.06 


- 


- 


38.9±7.0 


- 


051016B 


2.01±0.12 


- 


- 


0.047±0.012 


- 


060210 


2.25±0.05 


- 


- 


15.8±1.2 


- 


060604 


1.97±0.08 


- 


- 


1.31±0.19 


- 


060614 


2.01±0.09 


- 


- 


0.0010±0.0002 


- 


060908 


2.37±0.19 


- 


- 


26.9±7.1 


- 


060927 


2.09±0.16 


- 


- 


18.7±3.2 


- 
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Table 3: Spectral characteristics of the Swift BAT GRBs used to study the E. 
correlation. Fitting model is a single power-law N{E) oc E^^. The reference of redshift can 
be found at http : / /heasarc. gsfc.nasa.gov/docs/ swift/archive/ grhtoble/ . 



GRB 



BAT mean flux 
(15 - 150 keV) 



Reduced x^ 
(d.o.f) 



;i+^)^p 

(keV) 



E ■ 

-"^7,180 



ergj 



10 



^ erg cm ^ s ^ 



050319 


2.12±0.12 


Q 1O+0.63 


0.44 


;2i) 


- 


- 


3.24 


050814 


1.86±0.11 


o 94+0.15 


0.86 


;2i) 


712±146 


38.3±13.7 


5.3 


051016B 


2.55±0.22 


1 49+0.17 


1.13 


;i6) 


- 


- 


0.9364 


060202 


2.00±0.15 


-, n4+0.53 


1.51 


;i2) 


- 


- 


0.783 


060210 


1.67±0.05 


6.45±0.20 


0.73 


;2i) 


770±94 


43.8±9.6 


3.91 


060502 


1.32±0.04 


11.6±0.3 


1.19 


;56) 


339±53 


10.6±2.9 


1.51 


060526 


1.76±0.10 


3.99±0.26 


0.85 


;2i) 


480±90 


19.4±6.3 


3.221 


060604 


1.87±0.24 


q 90+0.50 

'J-^'J-0.49 


0.70 


;2i) 


396±173 


13.9±10.5 


2.68 


060605 


1.48±0.13 


0.818±0.067 


1.79 


;2i) 


199±63 


4.23±2.34 


3.78 


060607 


1.32±0.04 


6.75±0.15 


1.10 


;2i) 


635±96 


31.4±8.4 


3.082 


060714 


1.62±0.10 


2.88±0.18 


0.83 


;2i) 


304±71 


8.83±3.61 


2.711 


060729 


1.75±0.52 


0.609+°:?°^ 


0.60 


;i2) 


17.6 ±10.2 


0.065±0.065 


0.54 


060906 


2.26+°:^^ 


i.i4lS:i 


1.22 


;i2) 


- 


- 


3.686 


060908 


1.25±0.03 


13.4±0.3 


1.09 


;27) 


813±152 


44.9±11.3 


2.43 
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Table 4. Redshifts estimated by Lx — Tbrk relation. The reference of observed redshift can 

be found at http : / /heasarc. gsfc.nasa.gov / docs/ swift/ archive/ grhtoble/ . The pseudo 

redshift (*) obtained by another method can be found at 

http : / /cosmos.ast. obs-mip.fr/projet/catalogpZ.php. 



GRB Pseudo z (our results) 



Observed 



Pseudo z(*) 



050128 

050319 

050401 

050505 

050607 

050701 

050712 

050713A 

050713B 

050802 

050803 

050814 

050822 

050915A 

050922B 

051008 

051016B 

051109B 

060105 

060108 

060109 

060111B 

060115 

060203 

060204B 

060210 

060219 

060306 

060312 

060319 

060323 

060413 

060421 

060502 

060507 

060526 

060604 



-0.86 
-0.69 

+0.78 



0.49 
0.55 
0.46 



2.52; 

J^-O"-0.60 
^21+0.53 

2.12 
> 3.67 

^•''"''-5.36 
1 «7+l-22 
J^-"'-0.61 

9 CI +1.44 

^•"^-0.50 

+0.47 
-0.37 

1 ro+o.ao 

J^-^^-0.31 

+ 1.72 

-1.08 

1.29 

0.76 



1.97 



2 
2.49 



1 
1 

2.85 



qo+0.90 
•-'^-0.71 

QQ+0.54 

.OO_0.44 

+ 1.76 

1.09 

8.76;?;55g 

o co+3.82 
'J-^°-1.69 
7+6.19 
-3.41 



6.87: 



^■59tt:ll 



7.41 



+8.45 
-3.62 



+0.11 
0.11 



1.11 



5.46 



+ 11.53 
-5.46 
+7.39 
-2.79 



3.32 



+4.58 



1.81 

> 8.81 

Q 4-1+0.10 

^•^J^-0.11 

> 0.08 

1 01 +0-96 
^•^^-0.68 

7.871440 

Q ^0+3.36 
-^•"^-1.91 
-1.16 
-1.12 



3.24 
2.90 
4.27 

(z<5) 



(0.4-2.6) 

iz< 1.2, 1.71?) 

(5.3±0.3 (photometric)) 



0.9364 
(0.08?) 

(< 2.7) 



3.53 

(z<4) 
3.91 



3.83 



1.51 

3.221 
2.68 



z < 0.36, 5.2±2.2 



4.0±1.3 



3.1±1.1 



-28- 





Table 4- 


-Continued 




GRB 


Pscudo z (our results) 


Observed z 


Pseudo z(*) 


060605 


3.28+1-?? 


3.78 




060607 


0.681°:!^ 


3.082 




060708 


r or+6.92 
O.O0_2 81 


(z-1.8, z < 2.3) 




060714 


4.43l?:g 


2.711 




060719 


3.67l?:f, 






060729 


A cc+0.08 
^•^0-0.08 


0.54 




060805 


> 6.04 






060807 


o Qy+2.68 






060813 


1 77+0.38 

^■' '-0.33 




2.38±0.40 


060814 


i.93l^:S 


0.84 




060906 


4.691^2:1 


3.686 




060908 


7.9lll°2f 


2.43 




060927 


17.9Ii4.i 


5.47 




061021 


o 1 r+0.49 
^■-L<J_o.42 


(z < 2.0) 




061121 


1 /IO+0.30 
J^-4'3-0.26 


1.314 





